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en ' Abstract 
^ ' 

, The Ostrovsky equation is a model for gravity waves propagating down 

' a channel under the influence of Coriolis force. This equation is a 

modification of the famous Korteweg-de Vries equation and is also 
Hamiltonian. However the Ostrovsky equation is not integrable and 
, in this contribution we prove its nonintegrability. We also study local 

' bifurcations of its solitary waves. 
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1 Introduction 

^ • In this contribution we investigate the integrability and the local bifurcations 

^ of the solitary waves of the nonlinear equation 

{Ut - (3Uxxx + '2uux)x = lu, (1.1) 

where (3 and 7 are constant parameters. The equation Hl.l|l is known as the 
Ostrovsky equation |27j ! arid is often called the Rotation- Modified Korteweg- 
de Vries equation. It is a model for gravity waves propagating down a chan- 
nel under the influence of Coriolis force. In essence, u in the equation can be 
regarded as the fluid velocity in the x— direction. The physical parameter 

^On leave from the Institute for Nuclear Research and Nuclear Energy, Bulgarian 
Academy of Sciences, Sofia, Bulgaria. 
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7 > measures the effect of the Earth's rotation. The parameter (5 deter- 
mines the type of dispersion, namely /? < (negative-dispersion) for surface 
and internal waves in the ocean and surface waves in a shallow channel with 
an uneven bottom. The parameter value /? > (positive dispersion) applies 
to capillary waves on the surface of liquid or oblique magneto-acoustic waves 
in plasma. In the limit 7 = 0, (|1.H) apparently reduces to the integrable 
Korteweg-de Vries equation. More details about the Ostrovsky equation can 
be found in |31 IH", "^^ TU', fTT! . 

The equation has three known integrals of motion. Indeed, it can 
be formally rewritten in a conservation law form as follows {d = d/dx): 



'^t = (-{Ut- PUa:a:a: + 2uU:r)t) , (1-2) 

\7 / X 



2 4 

- P{2uul - T/^n,,) + p^u,{d-\) + y (a-2n)2 - y] ^. (1.4) 

From (|1.2j) - (|1.4|) we have the following three integrals, provided the solution 
u is in an appropriate class of functions, such that the integration (over the 
real line for Schwartz class functions, or over a period for periodic functions) 
makes sense (e.g. see for details): 

I[u] = j udx = 0, 7/0, (1.5) 

P[u] = ^ j u'^dx, (1.6) 

= /(fn2 + 2(0-V + y)dx. (1.7) 



Now we notice that the equation (ll.lj) can be written as 

_ d 6H 
dx 6u{x) ' 

where the symbol 5/5u denotes variational derivative. Moreover, (|1.<S|) can 
be further represented in a Hamiltonian form with a Hamiltonian H: 

ut = {u,H}. (1.9) 

The Poisson bracket is defined as 
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One can check that the bracket (|l.in|) is anti-symmetric and satisfies Jacobi 
identity. 

Also, we notice that the integral (|1.6)) has a meaning of a momentum. 
Indeed, it is related to the translation invariance of the Hamiltonian. Since 
H[u{x + e)] — H[u{x)] = 0, the expansion of / {H[u{x + e)] — H[u{x)])dx in 
e about e = gives (note that u{x) = 5P/5u) 



0= [JIL^dx= f ^ dx= [H P] - P, 

J 6u{x) dx J 6u{x) dx 6u{x) ' ' 

i.e. Pt = 0. The existence of (|1.5j) from the Hamiltonian viewpoint is related 
to the presence of the operator d/dx in the Poisson brackets. 

Apparently there is no second Hamiltonian formulation for the Ostrovsky 
equation, compatible with the presented one, i.e. this equation is not bi- 
Hamiltonian. Despite the existence of the three integrals, in Section 2 we 
demonstrate that the equation (|1.H) is not completely integrable for 7 7^ 0. 
In Section 3 we study the local bifurcations of its solitary waves. 



2 The Integrability Test 

In our analysis we use the integrability check developed in |^ |^ 0^ . 
(Another application of the method is presented in jl8j.) 

This perturbative method can be briefly outlined as follows. Consider 
the evolution partial differential equation 

ut = Fi [u] + F2 [u] + F3 [m] + . . . (2.1) 

where Fk\v\ is a homogeneous differential polynomial, i.e. a polynomial of 
the variables u, Ux, Uxx, d^u with complex constant coefficients, satisfying 
the condition 

Fk[Xu] = X''Fk[u], A G C. 

The linear part is Fi[u] = L(u), where L is a linear differential operator of 
order two or higher. The representation (|2.1j) can be put into correspondence 
to a symbolic expression of the form 

= nu;(6) + yai(ei,6) + ya2(6,6,e3) + -.. = (2.2) 

where a;(^i) is a polynomial of degree two or higher and ak{Ci,^,2, ■ ■ ■ S,k+i) 
are symmetric polynomials. Each of these polynomials is related to the 
Fourier image of the corresponding -Ffc[n] and can be obtained through 
a simple procedure, described e.g. in ^2^. Each differential monomial 
^no^rai _ _ (c}^n)"'9 is represented by a symbol 

,,ni//:0 tO el el e2 t2 pq \ 

\S1 • • • Sno^no+l ' ■ ■ ^no+ni^?io+»^i+l ' ' ' ^nQ+ni+n2 ' ' ' ^ml 
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where m = uq + ni + . . . + Ug and the brackets ( ) denote symmetrization 
over all arguments (i.e. symmetrization with respect to the group of 
permutations of m elements Sm)- 

(/(6>6> • • • )^n)) = ^ ^ f{S.a{l),Cai2), ■ ■ ■ ,Ca(n)) 

Also, for any function F (j2.2|) there exists a formal recursion operator 

A = 7] + uM^i,v) + u^MCuC2,r]) + . . . (2.3) 
where the coefficients </>m(Ci; ^2, ■ ■ ■(,m,'n) can be determined recursively: 



M^U^)=N^{Ci,v)^iaiiCi,v) (2.4a) 

(6 + 6 + • • • +^m)am(6,6,- • - ^^1,??) + 

m— 1 

^-^ \m — n + I 

n=l 

-(t>n{il, ■ ■ ■ in,ri)am-n{in+l, ■ ■ ■ im:V + il + ■■■ +Cn))} (2.4b) 

with 

?1=1 ?1=1 

and the symbols ( ) denote symmetrization with respect to Cii '^2, • • • Cm, 
(the symbol r] is not included in the symmetrization). 

We need also the notion of a local function, which can be defined as fol- 
lows. The function bm{^i,^2T ■ ■ im^'n) 1 ^ — 1 is called local if all coefficients 
bmn{^i,^2, ■ ■ ■ (.m), n = Ug, ris+i, . . . of its expansion as ^ 00 

00 

bm{^l,C2,---Cm,r]) = ^ 6mn(6,6, • • •Cm)^?"" (2.6) 

n=ns 

are symmetric polynomials. 

The integrability criterion can be summarized as follows |25j : 



Theorem 1. The complete integrability of the equation \2.1\) . i.e. the ex- 
istence of an infinite hierarchy of local symmetries or conservation laws, 
implies that all the coefficients \2.4\j of the formal recursion overator \2. S\) 
are local. 
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The equation can be written in the form 

ut = d'^ipd'^ + 7)u - 2uu^ (2.7) 

The symbohc representation of the operator is l/rj. Moreover, Theo- 
rem ^ can be apphed in this case as well [221 ■ The equation (|2.7|) can be 
represented in the form (|2.2|1 with 

w(Ci) = — (2.8a) 

«i(6,6) = -2(ei+6) (2.8b) 

Then from (fTH) : 

, ^ ^^HCl + V? q ^ 

"^'^^''"^^ 3/3ry2^f (6 + r?)2 - + r?)2 + ^^^i ^ ^ ^ 

</'2(6, 6, ^) = ^2,-3(6, 6)??'' + $2,-4(6, 6)^-^ + • • • (2.9b) 
The expansion of 0i (6 , rj) 1)2. 9a() with respect to tj is 

.)--A„-i_J^„-3 + J^„~4^27(l±Mi)„-5+ (2 10) 

and therefore there are obstacles to the integrability of since the co- 

efficients in this expansion are not polynomials in 6 for 7 7^ 0. 
For comparison, in the integrable case 7 = (KdV), we have 

M^I,V) = -^V-' (2.11a) 
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+^(ei + 6)(e? + ei)^-' + ... (2.11b) 

Therefore, the only completely integrable equation of the form (|l.lj) is 
the one with 7 = 0, i.e. the KdV equation. 



3 Solitary Waves and Local Bifurcations 

Solitary waves of the Ostrovsky equation of the form u{x, t) = (j){x — ct) = 
(j){z) satisfy the fourth-order ODE 

(l^zzzz - q(t>zz + p4> = , (3.1) 



where 



7 c 
z = x-ct, P=p' ^^'^"^ 
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Recently, it was proved by Liu and Varlamov |22j that solitary waves (j) exist 
if the speed c satisfies that c < 2-^/7^, /3 > with the zero mass, that is, 
j (j)dz = 0. Surprisingly, even if the mass of the solitary wave of the KdV 
equation is not zero, it is shown j2J that the limit of the solitary waves of 
the Ostrovsky equation tends to the solitary wave of the KdV equation as 
the rotation parameter 7 tends to zero. For /? < 0, solitary waves in the 
form of stationary localized pulses cannot exist at all [51 12UI I5T] . 

Equation (|3.1j) is invariant under the transformation z — > — z and thus 
it is a reversible system. In this section, we use the theory of reversible 
systems ^| El 1211 IMl 1^ 1^ to characterize the homoclinic orbits to the 
fixed point of ^'A.l^ . which correspond to pulses or solitary waves of the 
Ostrovsky equation in various regions of the (p, q) plane. 

The linearized system corresponding to (|3.H) 

(pzzzz - q(pzz +p(j) = (3.3) 

has a fixed point 

(j) = (t>z = (t>zz = (t'zzz = 0. (3.4) 

Solutions (f) = ke^^ satisfy the characteristic equation 

X'^ - + p = {), (3.5) 

from which one may deduce that the structure of the eigenvalues is distinct 
in four different regions of the (p, q) plane. These regions, as well as the 
eigenvalue structure, are shown in Fig.l. The regions are labeled as (1) - 
(4). The boundaries of these regions are the curves marked Co to C3 in Fig.l. 
In delineating the structure of the homoclinic orbits to fixed point (j3.4j) in 
various parts of the (p, q) space, we shall first consider the bounding curves 
Cq— C3 and their neighborhoods. Following this, we shall discuss the possible 
occurrence and multiplicities of homoclinic orbits to 1)3. 4() . corresponding to 
pulse solitary waves of the Ostrovsky equation in each of regions (1) through 
(4): 

a. Near Cq: This curve, on which the eigenvalues have the structure 
Ai_4 = 0, 0, ibA, A > and its vicinity have been considered in the context 
of reversible systems in jl6l I19j . In this region, a standard analysis yields 
the normal form on the center manifold 

Xi = X2, 

X2 = sign(^)Xi - ^^2, 

where fi is an unfolding parameter ^H]. For /i > 0, this yields a unique 
symmetric homoclinic solution 

Xi{t) = sech^ 
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in the vicinity of Cq. One may also show persistence of this homochnic 
solution in the original system ()3.4|) for fj, > |16j . 

b. Near Ci: Near Ci, which corresponds to the eigenvalue structure 
Ai_4 = 0, 0, ztzw, cj > 0, analysis of a four-dimensional normal from 
jl6j shows that on the side of Ci corresponding to Region 3 in Fig.l there 
is a sech^ homo clinic orbit. However, in Region 3, where the eigenvalue 
structure is that of a saddle-center Ai_4 = ±A, ±iuj, the fixed point H3.4() 
is non-hyperbolic. In fact it can be shown by the "antisoliton" method jHj 
that there are no soliton solutions of the system (|3.1|) for which the function 
4>{z) — > when |2:| — > oo together with its derivatives. 

c. Near C2: In this region, where Ai_4 = ziziuo, itio;, analysis of compli- 
cated normal form ^ EI shows the possible occurrence of so-called "enve- 
lope" homoclinic solutions of the form sech(A;t)e*"^ and with oscillating tails 
in the so-called "subcritical" form. However, occurrence or persistence of 
these solutions in the full nonlinear system (|3.1|1 is a non-trivial issue and 
each system must be analyzed separately [El El E] ■ 

d. Near C3: There is no small-amplitude bifurcation on C3, on which 
Ai_4 = ibA, ibA and the fixed point ()3.4|) remains hyperbolic. However, as 
we discuss below, there is a bifurcation across it causing the creation of an 
infinite multiplicity of homoclinic orbits. 

We turn next to each of the regions 1 to 4 in Fig.l to discuss the possible 
occurrence and multiplicity of homoclinic orbits in each region. 

Region 1. In this region Ai_4 = itA, ztiu; and the fixed point (|3.4j) is a 
saddle focus. Using a Shil'nikov type analysis, one may show jl3l I14j for 
general reversible systems such as (|3.1|) that the existence of one symmetric 
homoclinic orbit implies the existence of an infinity of others. Hence, we 
expect our system (|3.1|) to admit an infinity of such symmetric A^-pulses 
for each > 1. Here, a symmetric A'^-pulse oscillates N times in phase- 
space for z S (—00, 00) (or, more technically, crosses a transverse section to 
the primary symmetric 1 pulse times). In the context of the Ostrovsky 
equation, these would be A'^-peaked solitary waves, and we expect an infinite 
family for all > 1 for parameters (p, q) in Region 1 of Fig.l. 

Region 2. In this region Ai_4 = =bAi, ibA2 and the fixed point 1)3. 4|) is a 
hyperbolic saddle point. Thus, there is no a priori reason for multiplicity 
of homoclinic orbits in this region. However, depending on the actual form 
of the nonlinear term, a symmetric homoclinic orbit to 1)3. 4|) may exist (see 
[3 Ej ) . Also, depending on further conditions EHl , a further "orbit-fiip" 
bifurcation may cause complex dynamics in its neighborhood. In the context 
of our system (|3.1|) . these issues will need further investigation to establish 
possible existence of solitary wave solutions in this region of {p, q) space. 

Region 3. The generic situation in this region has already been considered 
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in the discussion above pertaining to the region near curve Ci. 

Region 4. In this region, Ai_4 = ziziiv, ^iuj2 and (|l-i.4|) is a focus. No homo- 
clinic orbits are known to exist in general here, although complex dynamics 
may occur jJl ISOl • 

Acknowledgments 

R. I. Ivanov and Y. Liu gratefully acknowledge the hospitality and support of 
the Mittag-Leffler Institute, Stockholm, where this research was performed 
during the semester program on "Wave Motion" in the Fall of 2005. R. I. 
Ivanov also acknowledges the funding from the Irish Research Council for 



Science, Engineering 


5 and Technology. 










J 
J 




< c f 







\ d") 








(4) 


<c \ 




(2) 













Co q 



(3) 

Figure 1 



8 



References 



Alexander JC, Grillakis MG, Jones CKRT and Sandstede B. Stability of pulses 
on optical fibers with phase-sensitive amplifiers. Z Angew Math Phys 1997; 
48: 175-192. 

Arnol'd VI and Sevryuk MB. Oscillations and bifurcations in reversible sys- 
tems, in: Nonlinear phenomena in plasma physics and hydrodynamics. Editor: 
Sagdocv R, Moscow: Mir, 1986. 

Benilov ES. On the surface waves in a shallow channel with an uneven bottom. 
Stud Appl Math 1992;87:1-14. 

Elphick C, Tirapegui E, Brachet ME, CouUet P and looss GA. A simple 
global characterization for normal forms of singular vector fields. Physica D 

1988;32:488. 

BuS^oni B, Champncys AR and Toland JF. Bifurcation and coalescence of 
a plethora of homoclinic orbits for a Hamiltonian system. J Dynam Diff Eq 
1996;8:221-279. 

Buffoni B and Groves MD. A multiplicity result for solitary gravity-capillary 
water waves via critical-point theory, preprint. School of Mathematical Sci- 
ence, University of Bath, 1996. 

Champneys AR and Toland JF. Bifurcation of a plethora of multi-modal ho- 
moclinic orbits for autonomous Hamiltonian systems. Nonlinearity 1993;6:665- 
721. 

Galkin VN and Stepanyants YA. On the existence of stationary solitary waves 
in a rotating field. J Appl Math Mech 1991 ;55: 939-943. 

Gilman OA, Grimshaw R and Stepanyants YA. Approximate and numerical 
solutions of the stationary Ostrovsky equation. Stud Appl Math 1995; 95:115- 
126. 

Grimshaw R. Evolution equations for weakly nonlinear long internal waves in 
a rotating fiuid. Stud Appl Math 1985;73:1-33. 

Grimshaw R, Ostrovsky L, Shrira V and Stepanyants Y. Long nonlin- 
ear surface and internal waves in a rotating ocean. Surveys in Geophysics 

1998;19:289-338. 

Gui G and Liu Y. Existence of the global low-regularity solutions for the 
Ostrovsky equation. J Funct Analysis, to appear 

Harterich J. Kaskaden homokliner Orbits in reversiblen dynamischen Syste- 
men. Master's thesis, 1993 

Harterich J. Cascades of reversible homoclinic orbits to a saddle- focus equi- 
librium. Physica D 1998;112:187-200. 

looss G and Kirchgassner K. Bifurcation d'ondes Solitaires en Presence d'une 
Faible Tension Superficielle. C R Acad Sci Paris Sr I Math 1990;311:265-268. 



9 



looss G and Kirchgassner K. Water waves for small surface tension: an ap- 
proach via normal form. Proc Roy Soc Edinburgh Sect A 1992;122:267-299. 

looss G and Peroueme M-c. Perturbed homoclinic solutions in reversible 1 : 1 
resonance vector fields. J Diff Eq 1993;102:62-88. 

Ivanov RI. On the integrability of a class of nonlinear dispersive wave 
equations. Journal of Nonlinear Mathematical Physics 2005;12:462-468; 
nlin.SI/0606046, 

Kirchgassner K. Nonlinearly resonant surface waves and homoclinic bifurca- 
tion. Advances in Applied Mechanics 1998;26:135-181. 

Leonov AI. The effect of earth rotation on the propagation of weak nonlinear 
surface and internal long oceanic waves. Ann N Y Acad Sci 1981;373:150-159. 

Levandosky S and Liu Y. Stability of solitary waves of the Ostrovsky equation, 
submitted. 

Liu Y and Varlamov V. Stability of solitary waves and weak rotation limit for 
the Ostrovsky equation. J Diff Eq 2004;203:159-183. 

Lombardi E. Homoclinic orbits to small periodic orbits for a class of reversible 
systems. Proc Roy Soc Edinburgh Sect A 1996;126:1035-54. 

Lombardi E. Orbits homoclinic to exponentially small periodic orbits for a 
class of reversible systems, Application to water waves. Arch Rational Mech 
Anal 1997;137:227-304. 

Mikhailov A and Novikov V. Perturbative symmetry approach. J Phys A 
2002;35:4775-90. 

Olver P and Jing Ping Wang. Classification of integrable one-component sys- 
tems on associative algebras. Proc London Math Soc 2000;81:566-586. 

Ostrovsky LA. Nonlinear internal waves in a rotating ocean. Okeanologia 
1978;18:181-191. 

Sanders J and Jing Ping Wang. On the integrability of homogenous scalar 
evolution equations. J Diff Eq 1998;147:410-434. 

Sandstede B, Jones CKRT and Alexander JC. Existence and stability 
of A^-pulses on optical fibers with phase-sensitive amplifiers. Physica D 
1997;106:167-206. 

Sevryuk MB. Reversible systems, Lecture Notes in Mathematics, 1211, New 
York: Springer, 1986. 

Varlamov V and Liu Y. Cauchy problem for the Ostrovsky equation. Discrete 
Dynam Systems 2004;10:731-751. 



10 



